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Three-dimensional stationary precession solitons with nonzero Hopf indices are found numericaUy 
by solving the LandauLifshitz equation. The structure and existence domain of the solitons are 
found. 
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Topological three-dimensional solitons attract considerable interest in many fields of physics including hydrody- 
namics, particle physics, cosmology, and condensed matter physics. In models with the three-component unit vector 
field n = ('^•1, '^•2, '^•3), where = 1, localized structures exist if the field n asymptotically approaches the vector no 
as |r| — > 00. Such fields map the U {00} space to the two-dimensional sphere and are classified by the homotopy 
classes 7r3(5'^) = Z and characterized by the Hopf invariant |Q] given by the expression: 
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where Fi — eyfen • (Vj^i x Vfcii) ^'^d curlA. = 2F. The invariant H admits the simple geometric interpretation as the 
linking number of two preimage closed curves corresponding to an arbitrary pair of points on the sphere. 

Stable three-dimensional solitons with H (the so-called knotted solitons) were studied numerically in the 
Faddeev-Niemi model H, ^ ||, |^, i.e., the nonlinear cr-model including terms with fourth-order derivatives. Defects 
with a nonzero Hopf invariant have been discussed in condensed matter physics since the pioneering works by Volovik 
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FIG. 1: Distribution of magnetization in a three-dimensional topological soliton with H — 3. The thin and thick curves are 
the preimages of the points (61, "l>i) — (O.Stt, 1) and (02, $2) = (0.27r, 3) on the sphere, respectively. 
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and Mineev for superfluid ^He Q and by Dzyaloshinskii and Ivanov |^ for a uniaxial ferromagnet. Simple reasoning 
based on the Derrick theorem || ] indicates the absence of nontrivial static threedimensional solitons with finite energy 
in the aforementioned media. However, dynamical structures of this kind, which are stabilized by the precession 
of the magnetization, can exist pO| , pA[ . Studies of threedimensional magnetic structures are of both academic and 
engineering interest in view of the development of new memory elements based on topological solitons in uniaxial 
ferromagnets p^ . 

Three-dimensional magnetic structures have been poorly studied until recently. A few original theoretical and 
experimental works on magnetic structures in superfluid ^He were published in recent years |l^, Spiral 

and cnoidal hedgehogs with H = were analytically described in in the Heisenberg model for an isotropic 

ferromagnet. Constant-velocity precession solitons with H = 1 were numerically found in the same model [ p^ . Both 
stationary (magnon drops) p9| and uniformly moving [^0| precession radially symmetric nontopological solitons in a 
uniaxial ferromagnet were analyzed. In addition, the moving topological soliton was recently analyzed pH |, but stable 
configurations with H were not found. 

In this work, we find dynamical structures, namely, stationary precession three-dimensional solitons with nonzero 
Hopf indices, in a uniaxial ferromagnet. Hereafter, the existence domain of such solitons is determined and their fine 
structure and main features are studied. 

The dynamics of the magnetization vector is described by the Landau-Lifshitz equation. In the case of negligible 
relaxation, this equation has the form 



dM _ 2^0 
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where the ferromagnet energy E = Eexch. + Eanis. is the sum of the exchange energy E^xch. and the magnetic 
anisotropy energy Eanis. with the parameter (3: 
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dr, (3) 



Eams. = ^ / + M/) dr. (4) 
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The energy of the magnctic-dipole interaction is negligible for large f3 values typical for many magnetic media. 
Using the parameterization 

rM = Afon, 

I n = {smQcos^ , sinQsin^, cosQ), 
we seek solutions describing stationary precession solitons of the form 

where ip is the polar angle of the cylindrical coordinate system (r, tp, z) and and 6 are the odd and even functions of 
respectively, and the z axis is the symmetry axis. It is known ^ that studying the dynamics of such structures 
is equivalent to solving a variational problem of minimizing the energy functional E for a fixed invariant of motion, 
N , equal to the total number of spin deviations: 

N ^ — / (1 - n^) dr = const. (7) 
2Aio J 

The precession frequency is determined by the formula 

The Derrick method applied to finite-energy field configurations yields the convenient formula for the precession 
frequency: 

1 E^xch. ^E^nis. /rt\ 
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In view of Eqs. (|), the energy functional takes the form 
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E = I / WEdrdz, (10) 
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Here, 



We = airMQ 



(11) 



where 

n* — {sin9cos(f), sin9sin(j), cosO). (12) 
The formula describing the Hopf index simpUfies drasticaUy for fields specified by Eqs. (^ ^ 

H = QT, (13) 

where the integer T is given by 
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To solve the variational problem numerically by the finite difference method, the energy functional was replaced 
by a function of the projections of the vector n* at the numerical-grid nodes. Since the field configuration for 
negative z values is determined by the symmetry of the problem, the discretization domain was the unit square 
0<r<l,0<z<l. 

The expected domain of soliton localization should be about the typical scale 

lo = J-s- (15) 



The parameters a and (3 were chosen to ensure that I/Zq ~ vD, where D is the number of grid nodes in one 
direction. The key point of this method is to ensure that the soliton is entirely contained in the simulation domain, 
but is sufficiently large to be adequately described by the discrete grid. In addition, to improve the performance, the 
numerical grid was nonuniform with the density of nodes increasing toward the coordinate origin. This means that 
the inhomogeneity was assumed to be localized near the center, while the field quickly reaches the asymptotic uniform 
state at infinity. The vector no = (0,0, 1) at the boundaries z = 1 or r — 1 of the domain corresponds to the ground 
state. 

The energy with an additive penalty function was minimized by the conjugate gradient method. A quadratic penalty 
function imposes a constraint on the norm of the vector n* which is unity at any point. Integral (^ is discretized to 
a linear form, and the imposed constraint is taken into account by applying a linear operator [p4[ in the calculation 
of the direction toward the minimal objective function for each iteration. 

A feature of our numerical analysis was the choice of the initial state, i.e., the trial field configuration from which 
the minimization begins. The vector corresponded to the ground state in the entire domain of simulation excluding 
a small quadrant near the coordinate origin. The vectors inside the quadrant lied in the orthogonal plane, and their 
directions varied randomly from point to point. Thus, the initial configuration is to a certain extent chaotic and does 
not even obey the criteria for the discretization of any smooth function. 

The result was directly tested in addition to the observed approach of the energy to its asymptotic value after 
many iterations. The frequency uj was calculated using Eq. The second derivatives and the discrepancy of 

Landau-Lifshitz equation (H) were numerically evaluated at each internal point of the simulation domain. 

For the solitons presented below, the numerical calculations based on Eq. ( p^ ) yielded T = 1. In this case, according 
to Eq. (p^, the Hopf index H = Q. Topological solitons of other types will be described elsewhere. The results of 
our present study are discussed below. 

Figure |l| shows a typical configuration of the field n for H — 3. Here, the coordinate axes are shifted for clarity and 
illustrate only the spatial orientation of the soliton. It is seen that the linking number of the preimages of two points 
of the S'^ sphere is equal to three. The toroidal surface corresponds to 9 = 7r/2. Note that n = no on the symmetry 
axis r = and n ^ no for r'^ + z"^ ^ oo. 
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FIG. 2: Contours of the angles parameterizing the unit vector n* for H — 3 and cu/uo = 0.5. The solid and dashed curves 
correspond to 9 — const and ctg{(j)) = const, respectively. 




FIG. 3: Distance dependence of the normalized energy density for solitons with the same precession frequency for various H 
values. The solid, dashed, and dashdotted curves correspond to H — 2, 3, and 4, respectively. 



The contours of the angles parameterizing the vector n* are plotted in Fig. ^ The curve = 7r/2 is a section of 
an axisymmetric surface of kind 1 (toroidal surface) in which a large fraction of the soliton energy is contained. The 
common point of the lines (f> = const is the south pole of the sphere S'^, at which = n. 

Analysis showed that it is most difficult to simulate solitons with small H values. A series of efficacious numerical 
results were obtained ior H — 3 and 4 using a 600x400 grid. However, in the case where H = 2, a, plausible result was 
obtained only for a 1600x800 grid. The cause is an interesting feature of the structure of the studied class of objects. 
This feature becomes apparent in numerical calculations of the normalized energy density 

1 f°° 

PE = WEdz. (16) 
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FIG. 4: Distance dependence of the normalized energy density of a soliton for several values of the precession frequency and 
H = 3. 



The dependences of p^; on r for various values of H are plotted in Fig. The maximum pE values are reached in 
the vicinity of 6* = tt/2. With a decrease in H, the localization domain corresponding to the first maximum of psir) 
becomes narrower and the maximum becomes steeper. The soliton structure becomes narrower, so that a finer grid 
with more nodes should be used. 



Three-dimensional solitons with H = exist for uj < O.OlScJo |19 , where ujq is the ferromagnetic-resonance 
frequency: 
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An increase in the precession frequency results in the compression of the soliton (Fig.H) implying the possible existence 
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of certain threshold ratio u) /ujq similar to the nontopological soliton jl^l • 

Figure || shows the dimensionless energy e = E/{ala) as a function oi uj/loq. It is seen that the soliton energy 
decreases with an increase in the precession frequency. A similar dependence for a magnon drop, calculated by the 
aforementioned procedure, is shown in this figure for a comparison. 

The discovered solitary structures with typical sizes of a few to a few tens of Iq are much smaller than cylindrical 
magnetic domains and, therefore, can potentially be used for information recording and storage if mechanisms for 
their generation and control will be developed in the same way as it is done now for magnetic vortices |l2|. This 
would provide the possibility of recording information in three-dimensional samples. 
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